arXiv:gr-qc/0503014v2 3 Nov 2007 


Reparameterization invariants for anisotropic Bianchi I cosmology 

with a massless scalar source 


D. C. Salisbury, J. Helpert, and A. Schmitz 
Department of Physics, Austin College, Sherman, Texas 75090-4440, USA 
E-mail: dsalisbury@austincollege. edu 
(Dated: November 2, 2007) 

Abstract 

Intrinsic time-dependent invariants are constructed for classical, flat, homogeneous, anisotropic 
cosmology with a massless scalar material source. Invariance under the time reparameterization- 
induced canonical symmetry group is displayed explicitly. 
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INTRODUCTION 


There is now general agreement in the relativity research community that time is a rela¬ 
tional concept, yet disagreements persist both on the context in which time is ’’recovered” 
and the role, if any, of the general covariance of Einstein’s classical theory of general rel¬ 
ativity. The contextual stances range all the way from the now largely dismissed claim 
that no true classical observable can vary with time (“frozen time”), to the expectation 
that variation in time will emerge only for appropriately selected quantum states from a 
timeless quantum theory of gravity. Attitudes regarding the underlying symmetry range 
from assertions that the full diffeomorphism symmetry group is lost forever in initial value 
formulations of general relativity, to our counterclaim that the full symmetry group is alive 
and well in canonical gravity. Indeed, our intention in this paper is to show how this canon¬ 
ical symmetry group can be used to introduce an observable relational time in a non-trivial 
gravitational model, an anisotropic universe with a massless scalar source. We will employ a 
canonically fixed intrinsic time, and we will compare and contrast our symmetry-based use 
of this time with several other approaches discussed in the literature. 

It is surprising that although it has long been advocated that intrinsically defined time 
may be employed in general relativity some researchers still believe that quantities that are 
invariant under the action of the full four-dimensional diffeomorphism group in canonical 
general relativity are constant in time, i.e., time is “frozen”. We believe this confusion 
persists in part because the full nature of the canonical symmetry group of general relativity 


has only recently been understood |l| . Crucial to this understanding is a distinction between 
gauge symmetry transformations and time evolution. Thus it is possible for the Hamiltonian 
to commute with invariants, which themselves commute with the individual constraints 
with which the gauge symmetry generator is constructed, and yet these invariants are not 
constants of the motion. 

The hrst proposal that material sources can be employed to intrinsically identify space- 
time events goes all the way back to Einstein as he first rejected and then ultimately em¬ 
braced the principle of general covariance in his general theory of relativity, as discussed in a 
classic paper by Stachel . The suggestion to use a spacetime scalar function of purely grav¬ 
itational dynamical variables as an intrinsic clock was hrst made by Geheniau and Debever 


, 14]- Komar and Bergmann showed subsequently that Weyl scalars could be expressed 


m 
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terms of canonical phase space variables and utdi^djo^x intrinsic coordinates in a classical 
Hamiltonian formulation of general relativity 


the use of intrinsic coordinates l9|. Rove 
in quantum gravity in the late 1980’s [10, 


li 

HI 


, 3, 8|. Dewitt was an early advocate for 

jegan to argue forcefully for their limited use 
12l |. Rovelli proposes a construction in which 


parameter(coordinate) and internal (intrinsic) time will in general not coincide. But the 
hxing of a parameter time in terms of a canonical intrinsic time is indeed the fundamental 
idea of this paper. Indeed, as we shall show, since our intrinsic time can be achieved through 
a time parameter transformation, every dynamical variable evaluated in the chosen special 
coordinate system is an invariant under time reparamterizations, and hence an observable. 

Marolf was the hrst to seriously apply the idea of intrinsic time to quantum cosmology 
. (See the review articles by Isham and Kuchar, and also a more recent paper by Smolin 
for the conventional analysis of the problem of time, including arguments for and against 
the generic existence of intrinsic time IJ, ll5|, Il6 |. 

It turns out that our cosmological model is exactly solvable with an arbitrary choice of 
time coordinate. We are thus able to demonstrate in detail that our phase space invari¬ 
ants do not change when the coordinate time is altered under the action of the canonical 
diffeomorphism-induced symmetry group presented by J. Pons, D.C. Salisbury and L. C. 
Shepley [l7|, [isl, ll9|, l20|. We construct invariants by choosing the evolving value of one 
of the cosmological expansion factors as an intrinsic time, i.e., by establishing correlations 
between this value and the scalar held plus the spacetime metric components, including 
the lapse function. As explained in detail in g the eelectton of a„ infnoeic ft.e cootdi- 
nate is equivalent to the establishment of a time-dependent gauge condition expressed in 
terms of phase space variables. Invariants under the action of the time-reparameterization- 
induced symmetry group discussed in the series of papers cited above can equivalently be 
constructed through the hnite action of this group along a gauge orbit. The lapse function 
in our cosmological model is one of the bases of this gauge transformation group. This is 
true since all coordinate time reparameterizations have a compulsory dependence on the 
lapse if we wish for them to faithfully reproduce in phase space the time-reparameterization 
symmetries present in conhguration-velocity space. The lapse must be intrinsically hxed, 
otherwise the gauge freedom is restricted as we shall show explicitly for the dynamical- 


variable-independent choice for the lapse made by Marolf 


13j. There are three essentially 


new features in our construction of invariants in cosmology. First, we carefully distinguish 
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between reparameterization symmetry (that requires the retention of the lapse as a dynam¬ 
ical variable) and coordinate time evolution. Second, we canonically generate inhnitesimal 
variations of our putative invariant phase space functions and show that these variations 
vanish. Along the way we will reveal in this non-trivial model how invariants can have 
vanishing Poisson brackets with the Hamiltonian, but can still have a non-trivial coordinate 
time dependence. And hnally, all of our invariants satisfy the original second order Einstein 
equations. 

In section 1 we give a brief overview of the general canonical symmetry group. We shall 
assume that coordinates have been chosen for our homogeneous cosmology such that all 
helds are independent of the spatial coordinates. Thus the only interesting diffeomorphism 
symmetry remaining is reparameterization of the coordinate time. We show in section 3 
that our coordinate choice may be made in the Einstein action, yielding a hnite dimensional 
Hamiltonian model. We hnd the general solution of the equations of motion in section 4. 
In section 5 we discuss the nature of the reparameterization-induced canonical symmetry 
group, and construct its generator. Then we construct variables that are invariant under the 
action of this generator. We conclude the section with a counting of the degrees of freedom 
of the model. 

The sourceless (Kasner) model has been quantized by Bojowald in the context of loop 
quantum gravity, as has the isotropic model with massless scalar source 0,Q]- Our analysis 
will be exclusively classical. Some loop quantum gravitational implications are addressed 
elsewhere. 


THE DIFFEOMORPHISM-INDUCED TRANSFORMATION SYMMETRY 
GROUP AND INTRINSIC COORDINATES 

We briefly summarize here the results of Pons, Salisbury and Shepley, and Pons and 
Salisbury, henceforth identihed respectively as (I) and (H). A group of transformations is 
said to constitute a symmetry if it leaves the form of the Lagrangian equations of motion 
unchanged. One consequence of this dehnition is that in order to be able to compare the 
form of the second order Einstein equations that follow from the Hamiltonian formulation of 
general relativity, we need to know how to transform both the dynamical variables and the 
coordinates under the action of the full diffeomorphism group. This means that we must be 
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continually cognizant of the relation between symmetry variations of the dynamical variables 
and the associated diffeomorphisms. Of course, another consequence of this dehnition is that 
under the symmetry transformation the full set of solutions to the second order equations 
of motion is transformed into itself. 

It has been standard practice, initiated by Dirac, to eliminate lapse and shift from the 
start as canonical variables. An immediate, obvious, and fateful outcome of this decision is 
that generically only the spatial diffeomorphism group survives as a symmetry. Obviously, 
with a fixed choice for time-dependent lapse and shift one is committed to a hxed time 
foliation of spacetime. And furthermore, it is apparent that the set of solutions (for this 
hxed choice of lapse and shift) does not include the complete set of solutions of the second- 
order equations of motion. The crucial observation here is that symmetries must transform 
the complete set of second order solutions into itself. If one insists in general relativity for 
example in keeping only the spatial metric and associated momentum as canonical variables, 
then one confronts the problem of how to construct a Hamiltonian in which the lapse and 
shift must vary (because we need to evolve the full set of solution trajectories), yet the lapse 
and shift cannot vary since they are not canonical variables. The unavoidable outcome of 
this analysis is that lapse and shift must be retained in the canonical formalism to have any 
hope of retaining the full diffeomorphism symmetry. 

But now we encounter another problem. If the lapse and shift are canonical variables 
we must insure that variations of these variables engendered by inhnitesimal coordinate 
transformations are faithfully reproduced as canonical transformations. This is the Legendre 
transformation projectability requirement. 

Consider, for example the variation of the lapse function N under an inhnitesimal coor¬ 
dinate transformation — e^. The resulting variation of the lapse is 

6N{t) := N'{t) - N{t) = iVe° + Ne^. (1) 

But since the Lagrangian does not depend on N it is possible to show that projectible 
functions in conhgur at ion/velocity space may not depend on N. Therefore e° must depend 
on N~^. Similar observation holds for time derivatives of the shift functions. It is then easy 
to show that projectable inhnitesimal spacetime coordinate transformations are of the form 

x'^ = x>^ - 6^C{x,g) ( 2 ) 
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where the arbitrary functions may depend explicitly on the coordinates x and non-locally 
on the spatial metric components and = {N~^, —N~^N°‘) is the normal to the hy¬ 
persurface of hxed coordinate time x^. It is important to observe that global translations 
in time, of the form x'^ = -|- constant are not projectable, i.e., time evolution in phase 

space cannot be put in correspondence with a Lagrangian symmetry. Since the projectable 
inhnitesimal diffeomorphisms have a compulsory dependence on the lapse N and shift N°‘ 
metric functions it is best to view the relevant symmetry group as a group of transformations 
on the metric and its time derivatives. Indeed, in the generic four-dimensional case, the full 
group depends non-locally on the spatial components of the metric due to the appearance 
of their spatial derivatives in the group Lie algebra. 

The Dirac Hamiltonian, the generator of time evolution, takes the form 

Hx = (3) 

where are arbitrary functions of spacetime coordinates. The P^ are the vanishing mo¬ 
menta conjugate to the lapse and shift, and the physical phase space is further constrained 
by the secondary constraints = 0. 

The projected inhnitesimal symmetry variations are generated in phase space by the 
generators 

G^it) = P^'^ + in, + , (4) 

where the structure functions are obtained from the closed Poisson bracket algebra 

{nA,nB'} = c^'b^Hc" , (5) 

and where spatial integrations at time t over corresponding repeated capital indices are 
assumed hereafter. The generators G^(t) act on phase space through the equal time Poisson 
brackets, and map solution trajectories into other solutions. In this sense, it is assumed 
that all phase space variables appearing in (jl]) are solution trajectories satisfying the Dirac 
Hamiltonian equations of motion. The generator G^{t) transforms solution trajectories into 
physically equivalent solution trajectories. Generally only one of these transformations in a 
given equivalence class is a time translation, namely the one for which happens to be equal 
to n^. We can generate all the members of an equivalence class by hxing a time foliation of 
spacetime, i.e., selecting initial lapse and shift functions. The action of the symmetry group 
then produces the full set of lapses and shifts, i.e., alters the time foliation. Thus although 
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it appears that the full diffeomorphism symmetry is lost in the canonical formalism, this is 
not the case. 

In fact, it was shown in (II) that it might be possible to transform physically inequiv¬ 
alent solutions to solutions that satisfy common coordinate conditions. In particular, it 
might be possible to hnd a global spacetime scalar function of the phase space variables 
that is a monotonically increasing function of the coordinate time for all physically distinct 
solutions. The same observation is made for spatial coordinates. Such global spacetime 
scalar functions are then chosen as intrinsic coordinates. All solutions in the intrinsic co¬ 
ordinate system are invariants under the action of the induced-diffeomorphism group, as 
we now demonstrate. Let y{x) represent phase space solution trajectories, and let S^{y{x)) 
represent four spacetime scalar functions, i.e., with the property that under the permissible 
coordinate transformations x'{x) 

S^{y'{x')) = S^{y{x)). (6) 

Dehne intrinsic coordinates through the coordinate transformation 

X^{x) ■.= S^{yix))=:a^ix), (7) 

with inverse 

x{X) = a-\X). (8) 

Next snppose that a coordinate transformation is nndertaken from x to x', and the transition 
to intrinsic coordinates is nndertaken from the primed coordinate system, i.e., 

X'^ix') := S'‘(!,'(x')), (9) 

Then it follows from ([6]) that the coordinate transformation from X^ to X'^ is the identity 
t r ansfor mat ion: 

X'>^{X) = S>^{y\x')) = S>^{y{x)) = a>^{x) = a>^{a-\X)) = X>^. (10) 

In other words, every geometrical object evalnated in the intrinsic coordinate system is 
invariant nnder coordinate transformations from x to x\ Indeed, this is not only a snfficient 
bnt also a necessary condition for invariance. 

We are able to constrnct invariants that exploit the diffeomorphism-indnced symmetry 
gronp. It can be shown in general that the following procednres are eqnivalent: 
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• Evaluate dynamical variables in the intrinsic coordinate system X^{x). 

• Impose the gauge conditions = S^{y{x)) and construct the corresponding Dirac 
bracket algebra. The idea here is that one and only one member of each equivalence 
class of solutions is selected with this gauge choice, namely the one for which the 
coordinate time and spatial dependence of the spacetime scalar functions are the trivial 
dependencies S^{y{x)) = x^. Generally each member will evolve under a different 
Hamiltonian, i.e., with different functions A^. However, all solutions will still obey 
the second-order Einstein equations. This gauge choice might not be globally, or even 
locally, implementable. It might eliminate entire equivalence classes. 

• Gauge transform phase space variables under the action of the diffeomorphism-induced 
group to solution trajectories that satisfy the gauge conditions above. The resulting 
gauge invariant variables satisfy the Dirac bracket algebra. 

We will illustrate all of these ideas with our cosmological model. 


CLASSICAL SYMMETRY REDUCTION AND HAMILTONIAN DYNAMICS 


We take the metric to be of the form 

/ 




n^it) 0 0 

0 a?'{t) 0 

0 0 lP‘{t) 


0 

0 

0 


( 11 ) 


y 0 0 0 c^{t) j 

where we have chosen to work in natural Planckian units, where ip, mp, and tp are respec¬ 
tively the Planck length, mass, and time. All of our variables, including the lapse variable 
n{t), are therefore dimensionless. We also take our space and time coordinates to be di¬ 
mensionless. We show first that the substitution of this metric into the general Einstein 
equations with massless source held yields the same equations of motion obtained 

by substituting this metric into the Einstein action and then varying the resulting reduced 
Lagrangian. 

The non-vanishing Ghristoffel symbols are 

h ha a 

T^U _ _ T^U _ _ _ _ 

■^00 ^1 ^ aa -*-00 „i 


n 




( 12 ) 



with similar expressions for a ^ b ^ c. (In the following we will not repeat this instruction). 
The resulting Ricci scalar is 


R = = 


2d 2b 2c 
a b c 


^ 2nd ^ 2nb ^ 2nc 2db 2dc 2bc\ 
na nb nc ab ac be J 

and the corresponding non-vanishing Einstein tensor components, where Ga /3 '■= Rap — 

2gctpR) 

^ db dc be _ 

Goo = -r + - + j-, 14 

ab ac be 

a% a?c a%n a^ch a%c 
““ bn"^ cn^ bn^ cn^ ben"^ 

Next we compute the non-vanishing components of the stress-energy tensor (which is 
obtained from the matter action — | / (Rx^/—g((),p) ■ We hnd 

Tm = r.„ = (16) 

Our Erst dynamical equation follows from the conservation law 

= 0 . ( 17 ) 

Only the a = 0 component is non-trivial, yielding 


Equivalently, 


dh (a b c 
n \ a b c 


The non-trivial Einstein’s equations G^^iy = are 

db dc be ■ 9 

-T + - + r = 2^*^ ’ 

ab ac be 


h (b c 
n\b c 


b c be -,9 

--— = 27r0 . 

b c be 


We next check to see whether we obtain the same equations of motion by hrst substituting 
our assumed metric into the Hilbert action and then varying this reduced action. The Hilbert 


action is 


d'^Xyf^R - - d^x^/^g^ 


— R'grav + 
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Referring to (fT^ and using = nabc the gravitational action becomes 


_ ^ 

'gray “ ^ 


S', = — I d'^x 


' d / 2hcd 2acb 2abc\ 2cdb 2bdc 2abc^ 
dt \ n n n j n n n ^ 


( 24 ) 


So ignoring the total time derivative, the effective gravitational action is 


Sgrav 


Also we hnd that 


h /■ 4 f dbc + abc + dbc' 

An J i n 


h f abccjP' 


Smatter I ^ ^ 

2 J n 


(25) 


(26) 


In the following we will ignore the spatial integral. (We could impose a spatial periodicity, 
in which case the Lagrangian would be multiplied by a dimensionless constant.) The total 
reduced Lagrangian becomes 


L = 


h / ; . . , . \ habc(fP‘ 


{dbc + abc + dbc 


+ 


(27) 


dvm ^ 2n 

Variation of this Lagrangian does indeed produce the equations of motion ([19]), fl20l) and 


Next we carry out the Hamiltonian analysis based on the reduced Lagrangian (j2711 . The 
canonical Hamiltonian is 


He = Paa + Pbb + PeC + P 4 ,^ + nn- L. 


The momenta conjugate to a, b, c, 0, and n are, respectively. 


dL 


h 


Pa = 


dd Aim 

dL habcc^ 


(he + 6c) , 


P<t> 


d<j) 


n 


and 


dL 

TT = — = 0. 

dn 


Solving fl29l) and fl30|l for the velocities, 

2im ( pb , apa p, 


a = 


c ^ be b J ' 


and substituting into fl28l) . we hnd 


n 


7rn ( 2paPb 2paPc 2phPc apl hpl cp. 


(28) 

(29) 

(30) 

(31) 

(32) 


+ -;^ + 


+ 


be ae ab 


+ 


npl 

2habe 


(33) 
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The Dirac Hamiltonian is 


Hx = H, + Att, (34) 

where A is a positive but otherwise arbitrary function of t. The resulting dynamical equations 


are 


dHx 27in ( pb , apa Pc 
a = —— = I-h 


dpa h \ c be b J ^ 

dHx TTTi (pI bpl cpl , 2pbPc\ , hnpl 

Pa - - j-1, n ol + 


da 


h \bc a?c a'^b 
■ dHx _ np^ 


a^bc ’ 


P<f> 


dpfj, Tiabc ’ 
dHx 


d(j) 
dHx 


= 0 , 


and 


TT = 


dn 


dHx ^ / 2paPb _ 2paPc _ 2pbPc ^ ^ 

dn h \ c b a be ac ab j 2habe 


Pt 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


Note that since tt is required to vanish (a primary constraint), preservation of this constraint 
over time yields a secondary constraint 

^2 ' 


H ■=- f ^ M ^ Q 

h \ e b a be ae ab I 2habe 


(41) 


GENERAL SOLUTION OF THE EQUATIONS OF MOTION 


The equations of motion ffT^ . (|20|) and can be easily solved by rewriting them in 
terms of the proper time r, dehned through the relation dr = n{t)dt. The general solution 
is 


and 


fl(t) — Oq 


b{t) = 6o 


c(t) = Co 


not + fo dt2 fo^ dtiA(ti) + tq 




To 


'^ot + fo dt2 fo^ dtiA(ti) + To 


To 

rt Jj. rt2 


not + fo dt 2 fo^ dtiA(ti) + tq 


To 


V 4,1 \ To 


(42) 

(43) 

(44) 

(45) 
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where tq is the value of r when t = 0. (Henceforth we let the variable symbol with subscript 
zero represent the value of the variable at time t = 0). The equations of motion ([19]), fl20|) 
and fl2ip result in the following conditions on the constants /cq, k/,, kc, and kf. 

ka + kb + kc= 1, (46) 


and 


h = ±\/l - iK + kl + kf). 

Notice that as a consequence of fITHD 

a{t)b{t)c{t) ao&oCo 


(47) 


(48) 


Thus we will be overcounting the complete set of solutions unless we £x this constant ratio. 
Let us therefore select 

To = ao^oCo, • (49) 

Let us also find the general solutions for the momenta. Substituting the configuration 
variable solutions flT3|) . (1441) . and time derivatives into fl29l) we hnd 

(50) 


h{l - ka) ( ao 

Po.[t) = -= PaO 


47ra(f) \a{t)) ’ 

where we recognize that the initial momentum pao is 

h(l - ka) 


PaO = 


Similarly, we hnd using fl5U]) and fl49]l that 


47ran 


(51) 


P<p{t) = 


\/4^ 


(52) 


REPARAMETERIZATION INVARIANTS 


Construction of reparameterization invariants 

The procedure presented in [l| is to hnd a function of canonical phase space variables 
that transforms as a scalar under the full diheomorphism group. In addition this scalar held 
must increase monotonically in parameter time for all physically distinct solutions in order 
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for it to define an intrinsic time. Of the many acceptable choices we will £x an intrinsic time 
T as, 

rp u\ ( ^ 0 ^ + -^0 ^^2 fo^ dtiA(ti) + aoboCo\ “ ^ ^ 

^ = ““ [ - 

This is to be nnderstood as a time reparameterization from time t to time T. Henceforth 
we will represent variables expressed in terms of intrinsic time, in addition to the intrinsic 
time coordinate itself, with capital letters. We constrnct invariants by carrying ont this 
transformation on our solutions. In the next section we will demonstrate explicitly that the 
resulting functions of phase space variables are invariant under the induced symmetry group 
corresponding to reparameterizations of the form t' = t — 

Let us hrst hnd the invariant expansion factor B{T) = ab{t{T)) where we recognized that 
6 is a scalar under the transformation fl53l) . It is not necessary to invert fl5^ for t{T) since 
the quantity appears in b{t) in the general solution flT3|) . Therefore 

we need only solve (1531) for this quantity: 


Therefore 


not + /o dt2 dtiA(ti) + apfepCo _ 

ao&oCo \aoJ 


/ rp \ kh/ka 

B{T) = b{t{T)) = bo(-) 
Vao/ 


(54) 


(55) 


Similar expressions are found for C{T). (Of course, A{T) = T.) Similarly, we hnd 



1 kj 


$(T)=0(t(T)) = 0o + y^^ln 

Also, since n(t) is a scalar density of weight one, 

n{t{T)) 


(56) 


N{T) = 


dr 

dt 


(57) 


Referring to fl55]) and fl5T)) . 


Therefore 


dT _ / ript + /o dt2 fo^ dtiA(ti) + apfepCo \ 

dt boCo \ o-oboCo j 

■ rp . l — l/ka 

n(t). 


boCo 

Oq 


ka-1 


n{t) 


kn 


boCo Vao 


N{T) = 


(58) 


(59) 
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Finally, since Pb(t) is a scalar, 


/an \ 

P,{T) = pb{t{T)) = pbo . 


(60) 


Summing up so far, we are to interpret the expressions (l55D . fl56|) . fl59|) . and fl60|) as func¬ 
tions of a time parameter T. The T dependence of our invariants is an explicit dependence, 
not an implicit dependence as an argument of the remaining dynamical variables. The in¬ 
variants depend only on the initial values of the remaining canonical phase space variables. 
Yet it turns out, as we will demonstrate shortly, that these initial values may be replaced 
by their evolved values at coordinate parameter time t. Indeed, any function A may be 
chosen in the Dirac Hamiltonian fIMl) . In other words, the dynamical-variable-dependent 
contribution to our invariants commutes with the Hamiltonian. Since the Hamiltonian con¬ 
tains the same constraints as the symmetry generator (j66j) . our putative invariants must 
indeed be invariants. Thus the fact that it is possible to replace the initial values by solution 
trajectories for arbitrary A is already a demonstration of invariance. Nevertheless, we will 
compute explicitly below the variations of our invariants generated by fl66l) . 

Indeed, substituting 


do — d(f) 


To 


not + fo dt2 fo^ dtiX{ti) -/ ao^oCo, 
and a similar expression for bo we hnd 


(61) 


B{T) = b(t) 



Similar substitutions into all of the invariant variables result in 


(62) 




N{T) = 


b{t)c{t) ( T 


kn 


and 


Pb{T)=pb{t) 


^a(t) ^ 

T 




^h/^a 


(63) 

(64) 

(65) 


These expressions will be used in section where we will demonstrate directly their invariance 
under the time reparameterization-induced symmetry group. 
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Finite transformations along gauge orbits 


Since we have at our disposal the full reparameterization-induced gauge symmetry group, 
we can alternatively use this group to construct gauge invariants. The idea is that given 
any solution of the equations of motion we can undertake a hnite gauge transformation to 
the unique solution satisfying the gauge condition expressing the fact that the intrinsic time 
has been chosen to be proportional to the scalar held. 

The generator G{^{t)) of inhnitesimal variations in phase space variables, induced by 
inhnitesimal reparameterizations t' = t — is according to ()T|) 


cm) 


^{t)H{t) + i{t)7r{t) 


/ 2pa{t)pb{t) _ 2pa{t)pc{t) _ 2pb{t)pc{t) 

h y c(t) b{t) a{t) 


plit) \ 

a(t)c{t) a{t)b{t) 2Tia{t)b{t)c{t) j 




a{t)pl{t) 

b{t)c{t) 


( 66 ) 


It is assumed in this expression that the phase space variables are solutions corresponding to 


the Hamiltonian fIMl) where a choice has been made for the function X{t). Given the hnite 
arbitrary function ^{t) it is easy to verify that the ehect of the hnite operator 


V^{s,t) = exp(^s{-, • (67) 

on phase space solutions yx{t) is to replace not + /g dG/o^ dtiA(ti) in the solutions corre¬ 
sponding to an initial choice of X{t) in the Hamiltonian by not -|- /g dt 2 /g^ dtiX{ti) + s^{t). 
Thus we obtain a one-parameter family of gauge transformed solutions of the form 


nsit) := V^{s,t)n{t) = n{t) + s^(f). 




0s (t) — 00 + 



not + Jo dt2 fo^ dtiX(ti) + tq + s^{t) \ 

To ) 

1 ^ f npt + /g dt 2 fo^ dtiX(ti) -|- Tg -|- s^(t) 


( 68 ) 

(69) 

(70) 


Our strategy is to hnd the descriptor ^(t, yx{t) that transforms a{t) to intrinsic time, and 
then to employ this descriptor in the hnite gauge transformation of all of the phase space 
variables. The resulting objects are manifestly invariant. Thus we set 


t = as=i{t) 

not /o dt 2 fo" dtiX(ti) rg -h ^(t; n(t), Og, bp, Cg) 




— Og 


To 


(71) 
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It is important to recognize that the left hand side of this expression is invariant under 
symmetry variations of all of the canonical variables. Changes in a{t), for example, corre¬ 
sponding to a given X(t), will be compensated by changes in ^(t-,n(t),ao,bQ,co). In other 
words, the time t will now be an independent parameter. 

Solving (ITTD for we hnd 


^(f;u(f),ao,6o,co) 



(72) 


According to fl68l) . fl69l) and (1701) the resulting gauge transformed invariants are precisely 
those obtained in fl59H . ([55]), and (l56|) . 

Demonstration of invariance 

Although in sections and we have already given rigorous indirect demonstrations 
of the gauge invariance of the phase space functions B(t), Pb(t), <I)(f)and N(t) under 
reparameterization-induced symmetry transformations, we will now check explicitly that 
they are indeed invariant under the inhnitesimal variations generated by (1661) for arbitrary 
inhnitesimal ^{t). 

We may take as our canonical variables y{t) at time t the solution pairs {a{t),pa{t)), 
{b(t), pb(t)), {c(t) , pc(t)) , and {(pit) , p^(t) = p^). Thus the variation of the solution trajectory 
b{t) at time t is 



(73) 


with analagous variations for b{t) and c{t). Referring to (l50|) and respecting the condition 


fH6|) we hnd that 



(74) 


Also, 


^Va{t) = {Pa{t),G{i{t))} 



a?{t)b{t) 27 ia‘^{t)b{t)c{t) 

H{t){ka - kl) 
Ana^{t)b{t)c{t) ’ 


hit)Pc{t) pl{t) _ b{t)pl{t) 
a^{t) b{t)c{t) a?{t)c{t) 


(75) 
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where in the last line we substituted the general solutions fl50D and took into account the 
conditions fH6l) and (HTj). In addition we hnd 




P<p 




ha{t)b{t)c{t) J 


(76) 


and 


Sp^ = 0. (77) 

Let us now check to see whether our phase space functions B(t), Pbit), $(t) and A^(t)are 
invariant under these inhnitesimal variations. Referring to ([62]), (TM]) . and (H6D we hnd 


6B{t) = {^b{t)a{t)-^^/^'^ - b{t)’^a{t)-^-’^»/^‘^6a{t)^ = 0. 

Similarly, referring to flM|) . flTill . and (H6|) . we hnd 

1 


SN{t) = —(a(f)i-^/^“5(6(f)c(f)) 

+ {b{t)c{t)) a(f)“^'^^“5a(t)^ = 0. 


Also, referring to (16^ . (ITil) . d75|) and (l50H . 


h 


SP,{t) = + !^a{t)-^+’^»/’^‘^pb{t)6a{t)j = 0 

And hnally, referring to fl63l) . fITill and fl76ll . 


(78) 


(79) 


(80) 


5$(t) = S(j){t) - = 0. 

/CaVdvr a{t) 


(81) 


Gauge fixing and counting of degrees of freedom 

As explained in detail in [l| the selection of an intrinsic time coordinate is equivalent 
to the choice of a time-dependent gauge condition. The gauge condition equivalent to our 
choice fl53l) is 

Xi = t — a{t) = 0. (82) 


The requirement that this condition be preserved under time evolution determines a second 
gauge condition; setting ^ = 0 we hnd 


0 = 


dXi 

dt 


dxi 

dt 


+ {Xi) -^a} 


k - 
1 

be 


(83) 
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Thus the second gauge condition is 


, . b(c)c(t) 

X2 = nit) - = 0. 


( 84 ) 


There are no more gauge conditions since the requirement that ^ = 0 merely hxes the 
function A. 

We investigate next the number of independent parameters at our disposal to vary the 
initial values of our phase space variables. Since there are five conhguration variables, a, b, 
c, n, and 0, the unconstrained phase space has ten dimensions. The two constraints tt = 0 
and (ITTD leave eight, and the two additional gauge conditions flH^ and flHTD should leave us 
six independent initial parameters in our invariant solutions. 

We note that all hve of the initial values Oq, bo, cq, uq, and 0o niay be freely specihed. 
However, according to (11^ and (jHTll . the only variable parameters in Pao, Pbo, and Pco are 
ka, kb, and kc- But these three parameters are subject to the conditions fHbll and fjTTl) . In 
addition, the momenta conjugate to $ and N are hxed by the other free parameters. Thus 
we do indeed hnd that the total number of independent initial values is six, corresponding 
to a dynamical system with three degrees of freedom. We may interpret the six independent 
parameters as labels of equivalence classes of solutions, where a single equivalence class is 
obtained by performing reparameterizations on solutions labeled by hxed values of the set. 

Our counting of degrees of freedom in this model appears to conhict with the counting 
irst undertaken by Ashtekar and Samuels in the Kasner model, in which there is no source 
24l. |25[ |. They pointed out that in the case in which the spatial topology is 3?^ the freedom 


to undertake global spatial diffeomorphisms reduces the phase space degrees of freedom 
to only one. We will address in a future publication the issue of homogeneity preserving 
spatial diffeomorphisms in the presence of sources. We observe that in any case our counting 
is correct when the spatial topology is that of a three-torus. Our focus in this paper has 
been on time reparameterizations. 


CONCLUSION 

As is the case with all generally covariant dynamical systems, the phase space version of 
the anisotropic cosmological model with a massless scalar material source is not covariant 
under changes in the time coordinate that depend only on the original coordinate time, i.e.. 
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coordinate transformations of the form 


t' = t\t). 


( 85 ) 


Nevertheless, there does exist a diffeomorphism-induced gauge symmetry that is realizeable 
as a canonical transformation group. The inhnitesimal elements of the diffeomorphism group 
are time reparameterizations of the form t' = t — where ^{t) is an arbitrary inhnitesimal 
function of the time parameter t, and n{t) is the lapse variable. The dependence on the 
lapse is required in order to produce variations of the metric and of the scalar held that are 
projectable under the Legendre map from conhguration-velocity space (the tangent bundle) 
to phase space (the cotangent bundle). The resulting phase space variations are canonical 
transformations, generated in the present model by fl66D . The phase space variables in this 
generator are solutions of the equations of motion corresponding to a choice of the function 
A(t), hence to an initial choice for the parameter time t. We should view t in this context 
as a parameter label identifying a one-parameter set of canonical variables. The canonical 
variables at time t generate on the same variables at time t a variation calculated via Poisson 
brackets of the same variables. The lapse solution n{t) is included among this one-parameter 
set of canonical variables. It must be included as a phase space variable in order to be able 
to implement the hnite diffeomorphism-induced symmetry group; successive inhnitesimal 
diffeomorphisms depend on the transformed lapse. 

We have employed this symmetry group in this paper to display in a non-trivial general 
relativistic model the time dependence of the metric and scalar variables that are invari¬ 
ant under diffeomorphism-induced symmetry transformations. The conventional argument 
against the existence of such variables is that if a variable cannot change its value under 
an arbitrary diffeomorphism t'{t), then obviously the variable must be constant in time. 
The problem with this argument is that t'{t) is not realizable as a phase space canonical 
symmetry group, as pointed out above. 

On the other hand, as noted in the Introduction, many authors have argued that non¬ 
trivial time evolution can occur when one employs some subset or combination of dynamical 
helds as a clock. One then seeks correlations between this dynamical clock and the remaining 


dynamical variables. Komar and Bergmann 


pioneered the use of intrinsic time in 


a canonical constrained Hamiltonian framework, stressing that intrinsic coordinates must 
be spacetime scalars. But they did not fully examine the canonical symmetry group. In 
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particular, they did not retain the lapse and shift as canonical variables. They were, however, 
the hrst to observe that an enlarged diffeomorphism symmet ry g roup existed, and that it 
possessed a compulsory dependence on the lapse and shift [401]. The resulting classical 
diffeomorphism-induced canonical transformation group has recently been studied, also in 
models in which additional gauge symmetries are present HQ,yQ- 

As we mentioned in the Introduction, many authors claim that the full diffeomorphism 
symmetry is lost in the transition to phase space. Dirac himself was led to “doubt how 
fundamental the four dimensional (symmetry) in physics is” and his remark stimulated Bar¬ 
bour’s engagement with timeless dynamical models. 2a, |27|] Kuchar has always maintained 
that the full diffeomorphism symmetry group is lost in conventional canonical gravity. We 
cite for example his remark in 1971 that the “canonical formalism necessarily destroys the 
spacetime covariance of the theory by cutting spacetime into slices” 2^ and a more recent 


observation with Kouletsis in an extension of a canonical covariant system that “Our exam¬ 
ination of the status of of spacetime diffeomorphisms in canonical description of covariant 
systems gives no support to a commonly held belief that spacetime diffeomorphisms are gen¬ 
erated by the Hamiltonian and momentum constraints”. 29(] Kuchar and his collaborators 
maintain that additional embedding variables must be included in order to recover general 


covariance in the canonical formulation of general relativity. 30, ImI, DO, l33|, [34 1 


Confusion over the role of the lapse is evident, for example, in Marolf’s treatment of 
time-reparameterization invariants in cosmology 13|. Although Marolf recognizes that the 
canonically implementable inhnitesimal time reparameterization symmetries are of the form 
t' = t — in posing the condition that variables transform as scalar densities of weight 
one under this transformation he implicitly requires that ^{t) is constant More re¬ 

cently, making the assertion that “translations in time are gauge symmetries”, Dittrich has 
employed a time independent lapse in her version of an implementation of the Rovelli con¬ 


struction of complete observables. |3g. 


37|, 


38l |. Since the lapse variable is generally not 


a gauge symmetry parameter, one is achieving a group average only in a very restricted 
context when integrating over it in the manner hrst advocated by Rovelli 4l|]. (Since spa¬ 
tial diffeomorphisms are realizable as canonical transformations without a dependence on 
the lapse and shift, it is possible that in contrast the shift functions might double as group 
parameters.) Rather, the descriptor ^{t) plays this role. In fact, in averaging over n{t) one 
is discarding signihcant physical information. We discuss in the Appendix the limited sense 
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in which time translation may be considered a gauge symmetry. 

As we have stressed throughout this work the full reparameterization covariance is re¬ 
tained only when the lapse is retained as a canonical variable and not viewed as an arbitrary 
coordinate function (and certainly not as a time-independent function). And the lapse in in¬ 
trinsic coordinates will generally have a non-vanishing Poisson bracket with other canonical 
variables. 

The classical, flat, homogeneous, anisotropic cosmological model with a massless scalar 
material source is exactly solvable with arbitrary time parametrization. We have chosen the 
value of one of the expansion factors as an intrinsic time. Referring to the general solutions 
of the equations of motion it is manifest that having fixed a set of six initial parameters as 
discussed in Section , for any value of the scalar field there corresponds unique values of 
the other phase space variables. The correlation is unique; it does not depend on the time 
coordinate (and correspondingly, the function A) with which one happens to be working. 
In other words, there is a unique map from the intrinsic time coordinate to the remaining 
field values. This means in particular that if we were to start with a coordinate time t 
and then transform to a new coordinate time t' before going to intrinsic coordinates, the 
transformation = T{t). In other words, the composite transformation from t to T' 

is the same as the original. Thus the symmetry transformation which transforms solutions in 
the coordinates t to solutions in coordinates t' can have no effect on the values of variables in 
the intrinsic coordinate system. This is the physical reason why our invariants, constructed 
as functions of phase space variables in the original time coordinate t, must be invariant 
under diffeomorphism-induced gauge transformations. 

Apart from our use of the full diffeomorpism symmetry group, we want to stress three 
essential differences between our method of construction of invariants, and the method em¬ 
ployed by Dittrich. 


36 


37( 1 First, every dynamical variable in our construction has an 


invariant form. The reason is that we are in fact merely transforming, via a dynamical-field- 
dependent coordinate transformation, to a special (intrinsic) coordinate system. Thus, even 
for just one choice of intrinsic coordinates, there are many more invariants than there are 
physically distinct solutions. The proof of invariance for the full set of dynamical variables 
requires that our intrinsic coordinate functions are true spacetime scalars. Kuchar and his 
collaborators have also identified problems that arise when the intrinsic time is a spatial 
scalar and not a spacetime scalar. 35| That our choice of intrinsic time is a spacetime 


21 







scalar can be demonstrated directly thongli the symmetry variation generated by fl66l) . cor¬ 
responding to the inhnitesimal time reparameterization t' = t — e{t), where e = ^/n. We 
hnd 


(5T(i) = = T(t)t(t), 


( 86 ) 


which is the transformation law of a spacetime scalar. 

Second, physically distinct solntions are identihed by their distinct initial conditions. 
For any choice of spacetime scalar T, there are many more Dirac observables (Dittrich’s 
terminology for invariants) than distinct initial conditions. 

Third, onr invariants will always satisfy the original second order Einstein eqnations. Of 
conrse, it will very likely not be possible generally to hnd global spacetime scalars with the 
correct coordinate dependence. It might perhaps be possible to patch fnnctions together. 
The “complete variables” considered by Rovelli and Dittrich will generally not satisfy the 
Einstein eqnations when the chosen “clock” partial variable is not a spacetime scalar. This 
raises deep conceptnal and even philosophical issnes. These two essential differences will 
featnre in a forthcoming extension of the intrinsic coordinate construction to the inhnite- 
dimensional cylindrical gravitational wave model. 

In an eventual quantum theory we must expect that the lapse will exist as an operator 
and that it will be subject to quantum huctuations, as was previously pointed out for the 
free relativistic particle Finally, we have stressed that evolution in the intrinsic time 
is non-trivial; time is not really “frozen” in general relativity. Though it may turn out that 
no suitable intrinsically dehned time exists in generic general relativity, where it does exist, 
as in this minisuperspace model, it is incumbent on us to explore its implications in an 
eventual quantum theory of gravity. This classical analysis s ugg ests a radical new approach 
to quantum cosmology, and to quantum gravity in general. |23|. 
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Appendix 


To the best of our knowledge this cosmological model has not been studied from the 
point of view of diffeomorphism invariants by either Rovelli or more recently, by Dittrich. 
However, they have clearly described their general program, and we shall implement it in 
this appendix so as to compare and contrast our results. 

In the language of Rovelli, our intrinsic time T is a “partial observable”. The objective is 
then to establish a correlation between the value of this intrinsically determined time and the 
values of the remaining dynamical variables. However, we encounter signihcant differences in 
the implementation of this program. First, neither Rovelli or Dittrich insist that their clock 
variable must always be a spacetime scalar. Second, we disagree in the implementation. 
Dittrich in particular considers exclusively evolution in proper time, corresponding to a 
gauge choice n = constant. This leads to inconsistencies, as we shall show. In addition, 
diffeomorphism invariants are then constructed out of constants of the motion. The resulting 
variables are manifestly invariant under global translation in time. That they are also 
invariant under the full diffeomorphims-induced group is also true, but this property is 
neither revealed nor explored by these authors. 

Dittrich identifies a phase space function T{y) as a clock partial variable. Letting f{y) 
represent any other independent phase space function, she constructs a corresponding one- 
parameter family of invariant phase space functions F'[/,r] ('r, y), where r is the parameter. 
The value of F'[/,t] for a hxed value of r is by dehnition the value that f{y) assumes when 
T{y) takes the value r. This formulation would be equivalent to ours at this stage if in 
addition T{y) were required to be a spacetime scalar. Indeed, Dittrich’s extension to held 
theory would also be equivalent to our introduction of intrinsic coordinates provided her 
additional “clock” variables were also spacetime scalars. They would then be interpreted as 
intrinsic spatial coordinates, in which case a better terminology might be “rod” variables. 
Henceforth we will assume that T{y) is a spacetime scalar, and we limit our remarks to the 
present cosmological model. 

We have shown that given any initial choice of time coordinate t it is possible to perform 
a canonically-generated symmetry transformation on the phase space function T{y{t)) such 
that, for all t, T{ys=i(t)) = t. In other words, for this presumably unique gauge r = t. 
Dittrich does not recognize the existence of the full gauge symmetry group. Rather, she 
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always works in a proper time gauge in which the lapse function n is a constant, so the 
arbitrary function A = 0 and the Hamiltonian is nif, where H is the secondary constraint 
fl34p ■ The only remaining gauge freedom in this case is global translation in the proper time, 
and it is generated by ^H, where is a constant. This is indeed a translation when acting 
on solutions in this special gauge since t' = t — ^ = t — constant. Furthermore, the functions 
F[f^T] are invariant under this gauge transformation. 

Let us now implement the Dittrich program. Inspecting the general solutions fl42|) - fl44|) 
we observe that a constant lapse results in a(t)b(t)c(t) = not, so let us choose T = abc. This 
gauge condition results in n = l.The resulting invariant functions in our notation and that 
of Dittrich are 

a{t) 


M'^d) = - 

{amt)c{t)) 

where we define /“(?/) := a. Similarly, 


■^Td = F^fo.^T=abc]i.'^D.yit)), 


(87) 


*h(r^) = 0(t) + \ —k^\n 

471- 


Td 


= Ft 




( 88 ) 


^a{t)b{t)c{t) 

Dittrich proves that her variables F are invariant under translations in t; for example, 

a{t + 5t) a{t) 

{a{t + 5t)b{t + 6t)c(t + {a{t)b{t)c{t))^‘" 

We have shown, on the other hand, that t is not restricted to proper time; solution trajec 


(89) 


tories for any choice of \{t) may used. And furthermore. 


a\{t) 


(ax{t)b^{t)cx{t)r 


is invariant under 


arbitrary diffeomorphism-induced symmetry transformations generated by 

We need to stress here that the intrinsic time uniquely fixes the Hamiltonian; both the 
lapse, and it’s time derivative, X{t) are determined by this choice. Therefore any choice 
other than T = abc will lead to inconsistencies in Dittrich’s program. To illustrate, let us 
pick r = a{t) and assume that t is the proper time. Then according to 




To ^ a{t) ^ ao (4) 


(90) 


Solving for A and substituting into the other proper time solutions we recover the invariant 


solutions 


and 


and 




kf,/ka 


In 


-) 

Oo/ 


(91) 

(92) 
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But there is one crucial difference. These expressions do not satisfy either the Hamil¬ 
tonian equations fl3^ - fl38l) . or the Einstein equations fl20|) - fl2ll) for n = 1. Rather, 
^ (^) ^ ” and \{td) = (^) ^ These inconsistencies do 

not arise in the intrinsic coordinate approach we have pursued in this work, where further¬ 
more it is recognized that the lapse is always an observable. The Dittrich procedure obscures 
this important point. Generally, the resulting lapse will have an explicit intrinsic time de¬ 
pendence, as well as a dependence on canonical variables, as is the case with fIMl) . The 
latter dependence will prove especially significant in quantum gravity since it will produce 
quantum fluctuations in proper time. 


[1] J. M. Pons and D. C. Salisbury, “The issue of time in generally covariant theories and the 
Komar-Bergmann approach to observables in general relativity,” Phys. Rev D71, 124012 
(2005) |gr-qc/0503013| 

[2] J. Stachel, “Einstein’s search for general covariance”, in Einstein and the History of General 
Relativity, D. Howard and J. Stachel, eds. (Birkhauser, Boston 1989) 

[3] J. Geheniau and R. Debever, Bull, classe sci. Acad. roy. Beige. 42, 114, 252, 313, 608 (1956). 

[4] J. Geh&iau and R. Debever, Helv. Phys. Acta Suppl. 4, 101 (1956). 

[5] A. Komar, Phys. Rev. Ill, 1182-1187 (1958). 

[6] P.G. Bergmann and A. Komar, “Poisson brackets between locally defined observables in gen¬ 
eral relativity,” Phys. Rev. Lett. 4, 432-433 (1960). 

[7] P.G. Bergmann, Phys. Rev. 124, 274-278 (1961). 

[8] P.G. Bergmann, “Observables in general relativity.” Reviews of Modern Physics 33, 510-514 
(1961). 

[9] B. DeWitt, “The quantization of geometry”, in Gravitation; an introduction to current re¬ 
search. L. Witten, ed. (Wiley, New York 1962). 

[10] G. Rovelli, “Is there an incompatability between the ways in which time is treated in gen¬ 
eral relativity and in standard quantum mechanics?”, in Conceptual Problems of Quantum 
Gravity, A Ashtekar and J. Stachel, eds. (Birkhauser, Boston 1991) 

[11] G. Rovelli, “Quantum mechanics without time: a model”, Phys. Rev. D42, 2638-2646 (1991) 

[12] G. Rovelli, “Time in quantum gravity: physics beyond the Schrodinger regime”, Phys. Rev. 


25 




D43, 442-456 (1991) 

[13] D. Marolf, “Quantum observables and recollapsing dynamics”, Class. Quantum Grav. 12, 
1199 (1995) |gr-qc/9404053 |. 

[14] C. J. Isham, “Canonical Quantum Gravity and the Problem of Time”, Presented at the 19th 
International Colloquium on Group Theoretical Methods in Physics, Salamanca, Spain, 1992 
gr-qc/9210011 

[15] K. Kuchar, “Time and the interpretations of quantum gravity”, in Proceedings of the 4th 
Canadian Conference on General Relativity and Relativistic Astrophysics: University of Win¬ 
nipeg, 16-18 May, 1991, G. Kunstatter, D.E. Vincent, and J.G. Williams, eds. (World Scien¬ 
tific, 1992) 

[16] L. Smolin, “The present moment in quantum cosmology: challenges to the arguments for the 
elimination of time”, |gr-qc/0104097 

[17] J. M. Pons, D. C. Salisbury and L. C. Shepley, “Gauge transformations in the Lagrangian 
and Hamiltonian formalisms of generally covariant theories,” Phys. Rev. D55, 658-668 (1997) 
gr-qc/9612037|. 

[18] J. M. Pons, D. C. Salisbury and L. G. Shepley, “Gauge transformations in Einstein-Yang-Mills 
theories,” J. Math. Phys. 41, 5557-5571 (2000) [gr-qc/99120861. 

[19] J. M. Pons, D. C. Salisbury and L. G. Shepley, “The gauge group in the real triad formulation 
of general relativity,” Gen. Rel. Grav. 32, 1727-1744 (2000) |gr-qc/9912087|. 

[20] J. M. Pons, D. G. Salisbury and L. C. Shepley, “Gauge group an reality conditions in Ashtekar’s 
formulation of general relativity,” Phys. Rev. D62, 064026-064040 (2000) |gr-qc/99120851. 

[21] M. Bojowald, ’’Homogeneous loop quantum cosmology” , Class.Quant.Grav. 20, 2595-2615 
(2003) [gr-qc/0303073|. 

[22] M. Bojowald, ’’Isotropic loop quantum cosmology”, Class.Quant.Grav. 19, 2717-2742 (2002) 
gr-qc/0202077|. 

[23] D. G. Salisbury, and A. Schmitz, “A generalized Schrodinger equation for loop quantum 
cosmology” |gr-qc/0702132 

[24] A. Ashtekar and J. Samuel, “Bianchi cosmology: the role of spatial topology”. Class. Quant. 
Grav. 8, 2191 (1991) 

[25] J. M. Pons, and L. G. Shepley, “Dimensional reduction and gauge group reduction in Bianchi- 
Type cosmology”, Phys. Rev. D58, 024001 (1998) |gr-qc/9805030 


26 


[26] P. A. M. Dirac, “The evolution of the physicist’s picture of nature”, Scientific American 208, 
45 (1963) 

[27] J. Barbour, The End of Time, (Oxford, 1999) 

[28] K. V. Kuchar, “Canonical quantization of cylindrical gravitational waves”, Phys. Rev. D4, 
955 (1971) 

[29] I. Kouletsis and K. V. Kuchar, “Diffeomorphisms as symplectomorphisms in history phase 
space: Bosonic string model”, Phys. Rev D65, 125026 (2002) 

[30] C.J. Isham and K.V. Kuchar, “Representations of space-time diffeomorphisms. 1. Canonical 
parametrized field theories”, Ann. Phys. 164, 288 (1985) 

[31] C.J. Isham and K.V. Kuchar, “Representations of space-time diffeomorphisms. 2. Canonical 
geometrodynamics”, Ann. Phys. 164, 316 (1985) 

[32] K. V. Kuchar and C. G. Torre, “Gaussian reference fluid and interpretation of quantum 
geometrodynamics”, Phys.Rev.D43, 419 (1991) 

[33] J.D. Brown and K.l V. Kuchar, “Dust as a standard of space and time in canonical quantum 
gravity”, Phys.Rev.D51, 5600 (1995) |gr-qc/9409001 

[34] J. Bicak and K. V. Kuchar, “Null dust in canonical gravity”, Phys.Rev.D56,4878 (1997) 
gr-qc/9704053 

[35] K.V. Kuchar, “Extrinsic curvature as a reference fluid in canonical gravity”, Phys.Rev.D45, 
4443 (1992) 

[36] B. Drittrich, “Partial and complete observables for Hamiltonian constrained systems” 
|gr-qc/0411013 

[37] B. Drittrich, “Partial and complete observables for canonical general relativity”. 
Class.Quant.Grav. 23, 6155 (2006) |gr-qc/0507106 

[38] C. Rovelli, “Partial observables”, Phys. Rev. D65, 124013 (2002) |gr-qc/0110035 

[39] C. Rovelli, Quantum Gravity, Cambridge University Press, 2004. 

[40] P.G. Bergmann and A. Komar, “The coordinate group symmetries of general relativity.” Int. 
Jour. Theor. Phys. 5, 15-28 (1972) 

[41] C. Rovelli, “The projector on physical states in loop quantum gravity”, Phys. Rev. D59, 
104015 (1999) [gr-qc/9806121| 

[42] J. M. Pons and D. C. Salisbury, “Invariants for cylindrical gravitational waves”, in preparation 

[43] D. C. Salisbury, ’’Gauge fixing and observables in general relativity”. Mod. Phys. Lett. A18, 


27 


2475 - 2482 (2003) |gr-qc/0310095| 

[44] Given that cj is a scalar density of weight one, it must transform under the canonically 
realizeable reparameterization t' = t — as 6 lo = ^ time derivative of ^ is 

missing in Marolf’s equation (3.10) 


